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Abstract 

Bright-rimmed clouds (BRCs) found in H II regions are probable sites of triggered 
star formation due to compression by ionization/shock fronts, and it is hypothesized 
that star formation proceeds from the exciting star(s) side outward of the Hll region 
(" small- scale sequential star formation"). In order to quantitatively testify this hy- 
pothesis we undertook BVI C photometry of four BRC aggregates. The amounts of 
interstellar extinction and reddening for each star have been estimated by using the 
JHK S photometry. Then we constructed reddening-corrected V/ V-I c color-magnitude 
diagrams, where the age of each star has been derived. All the stars turned out to 
be a few tenths to a few Myr old. Although the scatters are large and the numbers 
of the sample stars are small, we found a clear trend that the stars inside or in the 
immediate vicinity of the bright rim are younger than those outside it in all the four 
aggregates, confirming the hypothesis in question. 
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1. Introduction 

Bright-rimmed clouds (BRCs), or cometary globules, correspond to relatively dense 
clumps in a giant molecular cloud left unionized in the course of the expansion of the Hll 
region. So they are a sort of remnants of star formation activity in giant molecular clouds. 
But, at the same time, they are current sites of star formation. Theoretically, triggered star 
formation caused by the compression of the gas due to shock is expected to take place in 
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such a pre-existing dense clump in an H II region. This phenomenon is called radiation-driven 
implosion (RDI) and detailed model calculations were carried out by several authors (e.g., 
Bertoldi 1989, Lefloch, Lazareff 1994, Miao et al. 2006). Many molecular line and radio 
continuum observations of BRCs (Lefloch, Lazareff 1995, Lefloch, Lazareff, Castets 1997, De 
Vries, Narayanan, Snell 2002, Thompson et al. 2004a, Thompson et al. 2004b, Thompson, 
White 2004c, Morgan et al. 2004, Urquhart et al. 2006) show that the overall morphology 
and some physical properties of them are consistent with those of the RDI models. 

Actually, many BRCs are associated with signposts of recent and ongoing star formation 
such as Herbig-Haro objects and IRAS point sources of low temperature that meet the YSO 
criteria of, e.g., Carpenter, Heyer, Snell (2000). Sugitani, Fukui, Ogura (1991) and Sugitani, 
Ogura (1994) compiled catalogs (the so-called SFO Catalog) of altogether 90 BRCs associated 
with such IRAS point sources for the northern and southern skies, respectively. Subsequent 
near-IR (JHK) imaging observations of many of the BRCs in the above catalogs revealed that 
an elongated, small cluster or aggregate of YSOs are often associated with them (Sugitani, 
Tamura, Ogura 1995). Interestingly, "redder" stars tend to be located inside the BRC and 
closer to the IRAS source that corresponds to the reddest source and lies at the innermost end 
of the aggregate, and "bluer" stars tend to lie outside the cloud on the side of the O star(s) 
exciting the Hll region. These results led them to advocate the "small-scale sequential star 
formation (hereafter, S 4 F) hypothesis" that a BRC undergoes a few or several star formation 
events due to RDI always at its head part, which is located closer to the exciting stars(s) first, 
but then recedes outward in the Hll region as the the star formation proceeds, leaving groups 
of stars aligned in an age sequence. 

With the aim to strengthen the S 4 F hypothesis Ogura, Sugitani, Pickles (2002) made 
grism surveys of 24 BRCs for Ha emission stars. Altogether 460 Ha emission stars and 12 
Herbig-Haro objects have been detected in their vicinities. Presumably these Ha emission 
stars are mostly T Tauri stars (TTSs), although practically no follow-up observations have yet 
been done for them. These objects are found concentrated toward the head or just outside 
of BRCs on the side of the exciting star(s). One may suspect that there might be just as 
many such stars farther inside the BRCs, but concealed by the higher extinction. However this 
does not appear likely, because, as we see in Sect. 4, the Ay values of the Ha emission stars 
inside, e.g., BRC 14 are not much larger than those outside it. In addition in figures la to 
lz of Ogura, Sugitani, Pickles (2002) we do not find many faint (Ha emission) stars inside 
BRCs. So we think that the observed concentration of Ha emitters toward the heads of the 
BRCs supports the S A F hypothesis, although it may be exaggerated to some extent by the 
opacity effect. (Here we also add that as for the very rich AFGL 4029 IR cluster embedded 
in somewhat eastern part of BRC 14 (e.g., Deharveng et al. 1997) we believe that it formed 
spontaneously, not as a result of RDI; it is located too far away from the head of the BRC 14 
for the effects of RDI to propagate.) Recently Matsuyanagi et al. (2006) analyzed their deep 
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Table 1 . Log of observations 



BRC Date of Obs. Exposure Times (sec) 
x No. of Frames 



11NE 


05/9/27 


Ic- 


60x4 


V: 


120x4 


B: 


300x4 


12 


05/9/27 


ic- 


60x4 


V: 


100x4 


B: 


240x4 


14 


05/9/27 


le- 


50x4 


V: 


120x4 


B: 


300x4 


37 


05/9/25 


Ic- 


60x6 


V: 


200x6 


B: 


300x6 



Ji/if photometry of the BRC 14 region and showed that three indicators of star formation, i.e., 
the fraction of YSO candidates, the amount of extinction of all sources, and near-IR excesses 
of the YSO candidates, all showed a clear sequence from outside to inside the bright rim. This 
result further strengthens the S A F hypothesis. 

The best way to quantitatively testify the hypothesis is to estimate of the ages of the 
aggregate members and to compare them between different regions with respect to the bright 
rim. So we undertook BVI C photometry of BRC aggregates. The problems associated with 
this attempt, in contrast to the case of ordinary open clusters, are membership selection and 
reddening correction. We overcame these difficulties as described in Sec. 3, and constructed 
reddening-corrected color-magnitude diagrams (CMDs) for each BRC aggregate. The age of 
each star has been determined on the V / ( V-I c ) CMDs, and its spatial distribution has been 
examined. The result seems to confirm the S 4 F hypothesis. 

2. Observations and Data Reduction 

The imaging observations of the four BRCs in the BVI C bands were carried out on the 
2.0-m Himalayan Chandra Telescope (HCT) of the Indian Astronomical Observatory. HCT 
is located at Mt. Saraswati (altitude 4500m above see level) in the Himalayan region, but is 
remotely operated from CREST, Hosakote near Bangalore, via a satellite link. The instrument 
used is HFOSC in the imaging mode. The details of the site, HCT and HFOSC can be found 
at the HCT website 1 . The observations were made on 25 and 27 Sep, 2005 and a log of the 
observations is given in table 1. The sky conditions at the time of the observations were 
satisfactory with the seeing size of about 1.5". 

The data reduction was carried out at Aryabhatta Research Institute of Observational 
Sciences (ARIES), India. The initial processing of the data frames was made using the 
IRAF data reduction package. The photometric measurements were done using DAOPHOT 
II (Stetson 1987). The point spread function (psf) was obtained for each frame using several 
uncontaminated stars. The photometric accuracies depend on the brightness of the stars, and 

1 http : //www . crest . ernet . in/ iia/ iao/ iao .html 
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the typical DAOPHOT errors of our target stars ( V ~ 18 mag) are: 0.003 mag in I c , 0.005 mag 
in V, and 0.007 mag in B for BRCs 11NE, 12 and 37. However they are a little bit larger for 
BRC 14: 0.005, 0.005 and 0.007, respectively. The standard deviation of the residuals of the 
transformed V magnitude and B-V and V-I c colors of the standard stars were 0.010, 0.005 and 
0.012, respectively. Near the limiting magnitude of V ~ 20.5, the DAOPHOT errors increase 
to ~0.03 in all color bands. 

3. Star Selection and Reddening Correction 

BRC aggregates are very loose and composed of a small number of stars. In addition 
their Galactic latitudes are low, so there can be many field stars, although some of them may 
be associated with the Hll region to which the BRC belongs. Therefore we selected the targets 
from such BRCs that the associated aggregate contains relatively many stars. BRCs 11NE, 
12 and 14 belong to the Hll region IC 1848 and BRC 37 to IC 1396. For each BRC the Ha 
emission stars of Ogura, Sugitani, Pickles (2002) have been selected as the member stars. In 
the case of BRC 11NE and BRC 37, however, only a few Ha stars are available. But, as one 
sees in figures If and lp of Ogura, Sugitani, Pickles (2002), a relatively clear aggregate of 
stars can be noticed that appears to be associated with the BRCs; these aggregates show up 
more prominently in the 2MASS images. We, therefore, have added the stars having 2MASS 
counterparts from the aggregates, and they are indicated with "A" (additional stars) in the 
finding charts, figure 1. In BRCs 12 and 14 we deal only with the Ha stars of Ogura, Sugitani, 
Pickles (2002), so see the finding charts therein. 

Another difficulty is the extinction correction. The target stars must have different 
amounts of extinction caused by the dust in the BRC and/or Hll regions. Moreover, they 
are thought to be pre-main sequence stars and so have peculiar color indices, in particular, 
U-B ("UV excess"). Therefore the ordinary technique of U-B/B-V two-color diagram (2CD) 
to estimate the amount of reddening in open clusters cannot be used. B- V and V-I c are less 
peculiar in TTSs, but a new problem emerges: in the B-V/ V-I c 2CD the intrinsic line is nearly 
parallel to the reddening vector in the part for late-type stars, so this 2CD is practically not 
usable to estimate the amount of reddening. 

In order to overcome this difficulty we used the J-H/H-K s 2CD, where the Ay value for 
each star has been measured by tracing back to the intrinsic lines along the reddening vector 
found in Meyer, Calvet, Hillenbrand (1997). We adopted the intrinsic line of TTSs given by 
Meyer, Calvet, Hillenbrand (1997) and that of dwarfs of Bessell, Brett (1988). The former is 
used for Ha stars. The latter, used for non-emission stars, have been converted into the GIT 
JHK S system according to the equation given in Bessell, Brett (1988). For BRCs 11NE, 12 
and 37 the JHK S data have been taken from the 2MASS catalog and converted into the GIT 
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Fig. 1. Finding charts for stars in BRC 11NE and BRC 37. Stars with prefix "A" arc non-Ha emis- 
sion stars, and those without the prefix are Ha emitters of Ogura, Sugitani, Pickles (2002). They are 
reproduced from their figures If and lp, and are 2.5' by 2.5' wide each with north up and east to the left. 
In BRC 37 "jet" is part of the Hcrbig-Haro object HH 588 (see, Ogura, Sugitani, Pickles 2002) and the 
pair of thick tick marks indicate the nominal IRAS position. BRC 11NE does not harbor any IRAS point 
sources. 

JHK S system according to the relation given in the 2MASS web site 2 . Only stars with the 
photometric qualities 'AAA' have been adopted. For BRC 14 we have adopted the JHK S data 
used in Matsuyanagi et al. (2006), which were obtained with the infrared camera SIRIUS on 
the University of Hawaii 2.2-m telescope atop Mauna Kea. The identifications of the target 
stars between the HCT data and the 2MASS catalog or the SIRIUS data have been made on the 
basis of the J2000.0 coordinates within an accuracy of 2". From the Ay values thus obtained 
we have calculated the color excesses Ey_j c and Eb-v according to the ratio given by Munari, 
Carraro (1996). 

4. Results 

The photometric results for each star in each BRC aggregate are given in table 2. The 
second and third columns give the extinction-corrected V magnitude and V-I c color, respec- 
tively. The last column gives the location of the star with respect to the bright rim. Note that 
"outside/inside BR" (BR stands for bright rim) are the inside/outside of the Hll region. We 
cannot tell the accuracies of the values in Table 2, since the amounts of reddening corrections 
depend on several factors whose reliability is difficult to quantify. For each BRC a V /(V-I C ) 

2 http : //www . astro . caltech . edu/%7Ejmc/2mass/v3/transf ormations 
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CMD was constructed, where overlaid are the isochrones of Siess, Dufour, Forestini (2000) 
with Z = 0.02 and no overshooting, which are corrected for the distance given in table 2. 
These distances are taken from the SFO Catalog (Sugitani, Fukui, Ogura 1991). The age of 
each star has been determined according to these isochrones. In view of the limited page space 
we show as an example the Vq/(V-I c )o CMD for BRC 14 in figure 2. There are many stars 
which lie above the youngest isochrone of 1 Myr and their ages are estimated by extrapolation. 
The resultant ages appear to be reasonable, ranging from 0.1 to a few Myr (with only a few 
exceptions), which compare well with TTSs' lifetime of a few tenths to a few Myr. From this we 
would say that all the stars in our sample are probably aggregate members with practically no 
foreground/background stars mixed-in. The mass of the stars ranges from ~0.3 to ~2 M (for 
the 1 Myr isochrone), which is again reasonable. With these ages and masses the TTS nature 
of the Ha emission stars is almost confirmed. Vq/(B-V)q CMDs yield similar results. 

In each aggregate the locations of the stars with respect to the bright rim are grouped 
into two or three (only in the case of BRC 14) regions. The grouping is chosen so that there 
are at least four stars in each region with the exception of BRC 37 where we have altogether 
only six stars. The mean age and Ay of the stars in each region are summarized in table 3. It 
shows that the mean age of the stars inside or on the bright rim is always younger than that 
outside it, although the scatter is very large. The age difference is particularly significant in 
BRCs 12 and 14 where we have stars inside the bright rim. Moreover in BRC 14 the mean 
age constantly decreases from outside to inside the bright rim. Therefore our results appear to 
confirm the S A F hypothesis. 

From table 3 we see that stars inside the bright rim have larger Ay values than those 
outside it in the case of BRCs 12 and 14, although the scatters are large again. This confirms 
the result of Matsuyanagi et al. (2006) for BRC 14 (but with a smaller difference between the 
regions in our case). It is natural in view of the above general trend of the stellar ages. In the 
case of BRCs 11NE and 37, however, no difference is found between the two regions. 

5. Discussion and Conclusion 

Lee et al. (2005) argued that the presence of TTSs only near and/or outside the surface 
of BRCs would clearly indicate the evidence for triggered star formation, and in their study of 
BRCs in the Orion region they actually found TTSs only between the BRCs and the OB stars. 
In the present study we locate some TTSs inside the bright rim. But their Ay values show 
that they are not deeply embedded in the BRC. Therefore their apparent location may be due 
to a projection effect and they may actually be located near the front surface of the BRC, so 
in reality similar to the "on BR" stars. The difference between Lee et al. (2005) 's result and 
ours is probably due to the fact that our study is deeper and more small-scaled. 

The regions in each BRC aggregate show big scatters of the ages of the member stars in 
spite of a clear trend of the mean ages. The reason for this is not clear, but we suspect that it 
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Table 2. Results of dereddening and age estimation 



Star 


V 


(V-Jc)o 


(B-V) 


A v 


Age 


Location 










(mag) 


(Myr) 




BRC 11NE 


(1.9 kpc) 












Ha-6 


17.99 


1.84 


1.23 


1.46 


1.7 


on BR 


Al 


15.45 


1.57 


1.22 


1.41 


0.3 


" 


A2 


14.56 


1,57 


1.41 


0.00 


0.1 




A3 


17.44 


2.00 


1.40 


0.71 


0.7 


" 


A4 


16.25 


1.34 


0.80 


2.26 


1.5 




A5 


17.91 


1.75 


-0.21 


0.69 


1.9 


outside BR 


AG 


17.01 


1.67 


1.29 


0.78 


1.0 


" 


A7 


13.13 


0.85 


0.53 


0.00 


1.0 




A8 


19.13 


2.24 


1.65 


0.78 


1.9 




A9 


17.27 


1.61 


-0.65 


2.16 


1.8 




A10 


15.56 


1.58 


1.06 


3.09 


0.3 




All 


17.92 


1.8 


1.32 


1.55 


1.8 


)) 


A12 


16.84 


1.45 


1.11 


2.29 


1.9 




A13 


15.52 


1.54 


1.30 


0.00 


0.3 




A14 


16.26 


1.63 


1.21 


1.42 


0.6 




Ha-5 


18.32 


2.08 


1.52 


0.73 


1.4 




A15 


17.24 


1.82 


1.35 


1.38 


0.8 




Ha-4 


18.09 


2.05 


1.37 


0.61 


1.2 




BRC 12 


(1.9 kpc) 












Ha-1 


16.19 


1.61 


1.06 


2.61 


0.6 


outside BR 


Ha-5 


17.88 


1.52 


1.09 


0.00 


6.0 




Ha-10 


19.79 


2.22 


0.97 


0.00 


4.0 




Ha-17 


17.96 


1.88 


1.29 


1.06 


1.6 


)) 


Ha-19 


16.34 


1.70 


1.09 


0.94 


0.35 




Ha-09 


16.34 


1.95 


1.18 


1.62 


0.17 


on BR 


Ha- 11 


17.18 


2.10 




2.36 


0.35 


55 


Ha- 14 


14.81 


1.57 


0.99 


1.73 


0.1 


55 


Ha-16 


16.9 


1.83 


1.13 


1.41 


0.4 


)! 


Ha-20 


15.24 


1.29 


0.64 


3.54 


0.5 


inside BR 


BRC 14 


(1.9 kpc) 












Ha-5 


16.99 


1.84 


1.50 


2.18 


0.5 


outside BR 


Ha-12 


16.51 


1.46 


1.35 


3.40 


1.3 




Ha-16 


17.85 


1.85 




3.37 


1.5 




Ha-18 


19.07 


2.11 


7 


2.42 


2.5 




Ha-20 


17.29 


1.62 




3.14 


1.7 




Ha-23 


16.94 


1.39 




3.23 


3.0 




Ha-24 


18.93 


2.55 




2.14 


1.0 






12 3 

(V-Ic)o 



Fig. 2. Vq/ (V-I c )q CMD for the BRC 14 aggregate. Stars, which are all Ha emitters in this case, are 
plotted with different symbols according to their locations with respect to the bright rim. The isochroncs 
of Siess, Dufour, Forestini (2000) with Z = 0.02 and no overshooting are drawn with the correction for 
the distance of 1.9 kpc. The arrow indicates the reddening vector. 



Table 3. Summary of results 



BRC 


Region 


No. of 


Mean age 


Mean A v 






stars 


±st dev (Myr) 


±st dev (mag) 


11NE 


outside BR 


13 


1.26 ±0.60 


1.21 ±0.88 




on BR 


5 


0.86 ±0.71 


1.17 ±0.85 


12 


outside BR 


5 


2.51 ±2.43 


0.93 ±1.07 




inside & on BR 


5 


0.30 ±0.17 


2.13 ±0.86 


14 


outside BR 


7 


1.64 ±0.86 


2.84 ±0.57 




on BR 


5 


0.98 ±0.75 


3.39 ±0.66 




inside BR 


4 


0.50 ±0.22 


3.51 ±1.32 


37 


outside BR 


3 


1.93 ±2.22 


2.67 ±0.91 




on & near BR 


3 


0.92 ±0.63 


2.60 ±0.57 
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is caused by the proper motions of newly born stars. Jones, Herbig (1979) showed that young 
stars in the Taurus molecular cloud have velocity dispersions of 1-2 km/s in one coordinate. If 
this value is applicable to BRC aggregates, the stars must have moved 2 pc or so (in projection) 
from their birth place in 1 Myr. Since the size of the BRC aggregates is ~0.2 pc (BRC 37) 
to ~2 pc (BRC 12), this is sufficiently large for some of the stars which formed 1 Myr ago to 
move into the (projected) inside of the present bright rim. We suspect that many stars in H II 
regions have thier origin in already dispersed BRC aggregates except for the case where a rich 
cluster exists. 

The stars outside a BRC (i.e., inside the Hll region) have similar ages irrespective of 
the (projected) distance from the bright rim. It is our position that they essentially formed 
there inside the BRC that has since retreated to its present position. Then the above fact 
suggests that they were formed in the (first) collapse phase of the evolution of BRCs (Lefloch, 
Lazareff 1994). Presumably BRCs may have a second (or, even third) collapse phase due to 
their dumpiness or as the exciting O star(s) evolve(s) and increase(s) its (their) UV photon 
emission, or somehow they give birth to stars after this phase (i.e., in the cometary phase); 
and the stars on or near the bright rim, which generally are younger, may correspond to such 
origin. It is interesting to note that among the aggregate regions in table 3, those having the 
youngest mean age (i.e,, inside/on BR of BRCs 12 and 14) show the smallest age scatter. 

The age estimate from CMDs depends on the adopted distance; if we use a larger distance 
we get younger ages. For some BRCs in the SFO Catalog the distances are not well determined. 
However, since the results given in tables 2 and 3 are in a reasonable range and do not show 
any systematic differences between the BRCs, the distances adopted here are not in big error, 
presumably. Anyway, what is important in the present study is the relative ages among the 
stars in a BRC, and they are not changed by the error of the adopted distance. Another factor 
which affects the age estimation is duplicity of stars. On a CMD an unresolved binary star can 
be placed at an elevated position (by up to 0.75 mag), which results in a younger age. Some 
of the stars which fall above the 1 Myr isochrone (see figure 2) might be such stars. However, 
unresolved binaries must not be limited to this location, but exist among all the stars more or 
less equally on the CMDs. Therefore our main conclusion remain unchanged by this effect. We 
therefore conclude that our photometry of the BRC aggregates confirms the S A F hypothesis, 
although some problems remain. 

We thank the TAC and staff of HCT for the time allotment and for the support in the 
observations, respectively. Our thanks are due also to I. Matsuyanagi and the SIRIUS team 
for providing the JHK S data on BRC 14. We are grateful to the referee Dr. G.H. Herbig for 
his helpful comments. K.O. and A.K.P. acknowledge JSPS (Japan) and DST (India) for the 
finantial supports. Y.I. is supported by "The 21st Century COE Program: The Origin and 
Evolution of Planetary System" of MEXT. 
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